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The saturation model with DGLAP evolution is shown to give good description of the 
production of the charm and beauty quarks in deep inelastic scattering. The modifica- 
tions of saturation properties caused by the presence of heavy quarks are also discussed. 

1 Introduction 

The saturation model of Golec-Biernat and Wiisthoff (GBW) [2] has been very successful 
in describing both the inclusive, F%, and diffractive, F^ 3 \ structure functions of proton 
at low values of the Bjorken variable x. It incorporates the idea of parton saturation in a 
simple way by introducing the x-dependent saturation scale in the form Q 2 s {x) = (xq/x) x 
with the parameters xo and A determined from the fit to Fi . The model was then improved 
by Bartles, Golec-Biernat and Kowalski (BGK) [3] by including the DGLAP evolution of the 
gluon density, whose effects are important for the small-r part of the dipole cross section. 
This allowed to describe the new, more precise, HERA data. However, only the three light 
quarks contributions to the inclusive structure function F2 were considered by these authors. 

In order to consistently describe DIS one should take into account also the contributions 
of the heavy quarks since, as found by HI [5] and ZEUS [BJ, they may reach 30% for charm 
and 3% for beauty. In this short note Q] we present the results of the studies [I] of the 
DGLAP improved saturation model where also the charm and beauty quarks are present. 
The parameters of the model with the five flavors are fixed by the fit to the F2 experimental 
data. Then, the contributions of the heavy quarks to the proton structure function, F£ c and 
as well as the longitudinal structure function Fl and the diffractive structure function 
F^ 3 ^ are predicted. 

This study is related to that presented in [7] where the proton profile function is taken 
in a Gaussian form. In our case, however, the gluons are assumed to be evenly distributed 
over a certain area with a sharp boundary. Furthermore, in our approach we do not need to 
introduce an arbitrary mass for the light quarks as it was done in [3j [TJ 13 E] ■ 

2 The DGLAP improved saturation model 

The dipole picture of the photon-proton interaction at low x has been demonstrated to be 
a useful tool for calculating proton structure functions, both inclusive and diffractive. In 
this framework 7*p interaction is regarded as a two-stages process. Firstly, the photon with 
the virtuality Q 2 dissociates into a quark-antiquark par of a given flavor. This pair, referred 
to as a color dipole, is characterized by the fraction of the photon momentum carried by the 
quark (or antiquark) z and the qq transverse distance vector f. The splitting is described by 
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the photon wave function ^(r, z,Q , m^,e/), quantity fully calculable in quantum electro- 
dynamics (QED). In the second stage, described by the dipole cross section a{x, r), the color 
dipole interacts with the proton and, since the non-perturbative contributions are expected 
to be important, modeling of a{x,r) cannot be avoided. The expression for the inclusive 
proton structure function F2 may be, quite generally, written as 

F 2 (x,Q 2 ) = — I I dz\V f P (r,z,Q 2 ,m 2 f ,e f )\ 2 a(x,r), (1) 

47T Qf om j: p J JO 

where the sums run over the quark flavors / and the polarizations of the virtual photon P. 
In the BGK model the following form of the dipole cross section is used 

a(x,r) = <J jl -cxp ^--^—r 2 a s (^ 2 )xg(x,^ 2 )^ , (2) 

where fi 2 = C/r 2 + /ig. It interpolates between the GBW cross section [3] (at large r) and 
the perturbative result [TU] (at small r). Thus, both the feature of color transparency and 
gluon saturation are incorporated in this definition. The gluon distribution evolves with fi 2 
according to the leading order DGLAP equation, simplified further by neglecting quarks, 
with the MRST inspired initial condition 

xg(x,Ql) = A g x x *(l-x) 56 at Q„ = 1 GcV 2 . (3) 

Altogether, the model has five parameters ctq, C, /z 2 ,, A g and X g , which are determined by 
the fit to the F2 data. The fit with the charm and beauty contributions was performed using 
the recent data from HI [TT] and ZEUS [T2]. The HI data points were multiplied by 1.05 to 
account for slightly different normalizations between the two experiments. Since the dipole 
picture is expected to be valid only at sufficiently small x we restricted ourselves to x < 0.01. 
Due to the fact that the gluon density is evolved according to DGLAP equations the model 
is supposed to work for any value of photon virtuality. Thus, in the fit, we used the entire 
range of Q 2 covered by the data. This gave us 288 experimental points. The light quarks 
were taken massless and the typical values for the masses of the heavy quarks were used, 
namely m c = 1.3 GeV and mb = 5.0 GcV. The number of active flavors was set to 4 (for 
the fit with charm) or 5 (for the fit with charm and beauty), the value of Aqcd = 300 MeV, 

and the argument in the dipole cross section was modified x — ► x (l + 4m 2 /Q 2 J similarly 

to 0[3]. 



3 Fit results and model predictions 

The results of the fit with heavy quarks are summarized in Table [TJ where also the light 
quarks fit parameters from [3] are recalled for reference. We see that the quality of the fit 
in terms of % 2 /ndf is good. Adding heavy quarks results in a rather dramatic change of the 
parameters of the model. In particular, the sign of the power X g alters which means that the 
initial gluon distribution grows with decreasing x oppose to the the case of the light quarks 
fit where it was valencelike. 

The predictions for the heavy quark contributions to the structure function, F 2 CC and 
are presented in Fig. [TJ We observe very good agreement with the data points from 
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X^/ndf 


light + c + b 


22.7 


1.23 


- 0.080 


0.35 


1.60 


1.16 


light + c 


22.4 


1.35 


- 0.079 


0.38 


1.73 


1.06 


light 


23.8 


13.71 


0.41 


11.10 


1.00 


0.97 



Table 1 : The results of our fits with heavy quarks to the F2 data and the massless fit from [3] . 




Figure 1: Predictions for the charm (x4™) and beauty (x8 n ) structure functions from the 
BGK model with heavy quarks compared with the HI data [5]. 



HI [5]. This persists even for x > 0.01 that is beyond the range used in the fit to determine 
the model parameters. 

The diffractive structure function was also calculated and good agreement with 

the HI and ZEUS data was found. Likewise, the longitudinal structure function Fl obtained 
from our analysis agrees with the HI estimations. For more details on Fl and F® the 
reader is referred to [I] . 

Let us finally discuss the effect of heavy quarks on the position of the critical line. This 
line in (x, Q 2 ) plane which marks the transition to the saturation region is plotted in Fig. [21 
We have checked [1] that the presence of heavy quarks shifts this line slightly towards the 
smaller values of Q 2 at low x. Similar behavior has been already observed in the GBW 
model [2] . Let us also point out that the critical line presented in Fig. [2] is very similar 
to that obtained by Soyez [9] in the modified Iancu, Itakura and Munier (IMM) model [8j 
with heavy quarks (see [4] and [9] for the precise, slightly different, definitions of the critical 
lines). 
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4 Conclusions 



We have shown that the DGLAP improved sat- 
uration model with heavy quarks provides suc- 
cessful description of a considerable number of 
quantities measured in DIS. In particular, it pre- 
dicts correctly the charm and beauty contribu- 
tions to the proton structure function. This re- 
sult is quite remarkable given the simplicity of 
the framework we use. This may suggest that 
the kr factorization is a more efficient way of 
approaching DIS at small x (see also [13]) or be 
considered as an argument supporting the idea 
of saturation at HERA. 
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